Astronomy & Astrophysics manuscript no. rsclucamshells05 


©ESO 2010 


March 1,2010 





High-resolution HST/ACS images of detached shells around 

carbon stars 

H. Olofsson 1 ' 2 , M. Maercker 2 , K. Eriksson 3 , B. Gustafsson 3 , and F. Schoier 1 



o 

(N 



in 

6 

5— ( 



> 

(N 

m 
q 

cn 
o 
o 



X 



1 Onsala Space Observatory, Dept. of Radio and Space Science, Chalmers University of Technology, SE-43992 Onsala, Sweden 
e-mail: harts . olof sson@chalmers . se 

2 Department of Astronomy, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden 

3 Department of Physics and Astronomy, Uppsala University, SE-75 120 Uppsala, Sweden 



Received 21 December 2009; accepted 28 January 2010 



ABSTRACT 



Context. Overall spherically symmetric, geometrically thin gas and dust shells have been found around a handful of asymptotic giant 
branch (AGB) carbon stars. Their dynamical ages lie in the range of 10 3 to 10 4 years. A tentative explanation for their existence is 
that they have formed as a consequence of mass-loss-rate modulations during a He-shell flash. 

Aims. The detached shells carry information on their formation process, as well as on the small-scale structure of the circumstellar 
medium around AGB stars due to the absence of significant line-of-sight confusion. 

Methods. The youngest detached shells, those around the carbon stars R Scl and U Cam, are studied here in great detail in scattered 
stellar light with the Advanced Survey Camera on the Hubble Space Telescope. Quantitative results are derived assuming optically 
thin dust scattering. 

Results. The detached dust shells around R Scl and U Cam are found to be consistent with an overall spherical symmetry. They have 
radii of 1972 (corresponding to a linear size of 8xl0 16 cm) and 777 (5xl0 16 cm), widths of 172 (5xl0 15 cm) and 076 (4xl0 15 cm), and 
dust masses of 3xl0~ 6 and 3xlO~ 7 M Q , respectively. The dynamical ages of the R Scl and U Cam shells are estimated to be 1700 and 
700 yr, respectively, and the shell widths correspond to time scales of 100 and 50 yr, respectively. Small-scale structure in the form of 
less than arcsec-sized clumps is clearly seen in the images of the R Scl shell. Average clump dust masses are estimated to be about 
2xl0~ 9 M . Comparisons with CO line interferometer data show that the dust and gas shells coincide spatially, within the errors (< 1" 
for U Cam and » 2" for R Scl). 

Conclusions. The results are consistent with the interpretation of geometrically thin gas and dust shells formed by a mass-loss eruption 
during a He-shell flash, and where interaction with a previous wind plays a role as well. The mass loss responsible for the shells must 
have been remarkably isotropic, and, if wind interaction plays a role, this also applies to the mass loss prior to the eruption. Clumpy 
structure is present in the R Scl shell, possibly as a consequence of the mass loss itself, but more likely as a consequence of instabilities 
in the expanding shell. 

Key words. Stars: AGB and post-AGB - Stars: carbon - Stars: evolution - Stars: late-type - Stars: mass loss 



1. Introduction 

Mass loss from the surface is an important characteristic of stel- 
lar evolution on the asymptotic giant branch (AGB). It is a com- 
mon property of most M-type, S-type, an d all C-type AGB stars 
that has been well esta blished for decades dOlofsson et aU[l993l 
iRamstedt et alll2009l) . Yet, many of its details, such as the mech- 
anism behind it, and its evolution with time and e.g. dependence 
on stellar mass, are essentially unknown and remain the major 
obstacle for understanding stellar evolution on (and beyond) the 
AGB in detail as well as the cont ribution that AGB stars make to 
the galactic chemical evolut ion (iHabingl Il996t IWillsonl l2000t 
ISchroder & Sedlmavril200lt) . 

It is particularly important to understand the temporal evo- 
lution and the dependence on direction of the stellar mass loss. 
The former determines to a large extent how the star evolves, 
while the latter has a profound effect on the circumstellar evo- 
lution beyond the AGB, e.g., the formati on of planetary neb- 
ulae. Geometrically thin detached shells dSchoier et al.L l2005h 
are a phenomenon that bears on both issues. Some carbon stars, 
of which less than ten are known, show this phenomenon. It 



has been suggested that these shells are the result of strong 
mass-loss modulations during a thermal pulse ( Olofsson et all 
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ll990HSchroeder et allll998l:lwachter etaUl2002|) and that fhev 
are additionally affected by an interacti on with the surrounding 
(relic) circumstellar enve lope (CS E) (|Steffen & Schonberneti 
l2000HSchoier et al.Ll2005h . Recentlv. lMattsson et all d2007[) pre- 
sented models where the response to structure variations dur- 
ing a thermal pulse of the dynamical atmosphere and the ex- 
panding gas and dust are studied in some considerable detail. 
Geometrically thin shells appear under certain circumstances as 
an effect of a mass-loss eruption and a subsequent interaction 
with a previous slower wind. Remarkably, these shells imply 
that at le ast during this phase, the mass loss is very close to 
isotropic dOlofsson et all 1 19961 120001: iGonzalez Delgado et all 
120031: iMaerckeret all 120101) . The study of thin, detached shells 
also has a bearing on the small-scale structure of the circumstel- 
lar medium, since the line-of-sight confusi on is limited by the 
thinness of the shells dOlofsson et alll2000h . 

Most of the information on these detached shells stems 
from CO radio line obs e rvatio ns. The first detections were 
made by lOlofsson et al.l d!988l) . and to this date there are 
seven known carbon stars with this type of shell dSchoier et all 
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120051) . iLindqvist et all (1 19991) and lOlofsson et all (l2000t) used 
the IRAM PdB mm-wave interferometer to show that the 
shells are geometrically thin (width/radius < 0.1) and remark- 
ably spherically symmetric. ISchoier et alJ d2005l) found evi- 
dence of effects of interacting winds, i.e., the shells are af- 
fected by their progress in a previous slower stellar wind. 
IGonzalez Delgado et al.l (1200 ll) showed for the first time that 
these shells could be ima ged in stellar light scattered in the 
circumstellar medium , and IGonzalez Delgado et all ( 120031) and 
iMaercker et alJ d2010h . following up on this using polarimet- 
ric imaging, showed that both dust and atoms act as scat- 
tering agents. No M-type or S-type AGB star was found 
with sim ilar geomtrically thin shells, despite extensive s earches 
(see e.g jKerschbaum & Qlofssonlll999t(Ramstedt et al.Ll2009t) . 
Detached dust shells around AGB and post-AGB objects have 
been seen as well, but with much coarser r esoluti o n ob- 
servations .Waters etall Il994t llzumiura et alJ, 1 19961 Il997t 
lHashimoto et al.Lll998tlSpeck et al.Ll2000t) . 

We note here that some detached shells around AGB stars 
may have a different origin. Neutral hydrogen 21cm observa- 
tions of extended CSEs show that some AGB stars are sur- 
rounded by large detached shells, whose presence is most 
likely due to interactions between circumstellar winds of dif- 
ferent epoch s or an interaction with the surrounding interstel- 
lar medium (iLibert et all I2007J). Similar conclusion s are drawn 
based on Spitzer observations dWareing et al., 200a). 

In the CO mm observations bv lOlofsson elaTl (1 19961 l2000T) 
it was found that the detached shells showed a clumpy struc- 
ture. The optical observations bv lGonzalez Delgado et alJd2001l) 
verified this finding optically to some extent, but the angular 
resolution was severely limited by seeing. Moreover, scatter- 
ing in the terrestrial atmosphere limited the studies of the shells 
close to the star. In order to increase the resolution and study 
the circumstellar envelope closer to the star, we decided to use 
the Hubble Space Telescope (HST). We present broadband fil- 
ter images of the circumstellar environments of the carbon stars 
R Sculptoris (R Scl) and U Camelopardalis (U Cam) obtained 
with the Advanced Camera for Surveys (ACS) on the Hubble 
Space Telescope (HST). 



2. Observations 

The carbon stars R Scl and U Cam (their stellar properties are 
summarized in TableQ]) were observed with the ACS on the HST 
in November 2004 and January 2005, respectively. The images 
were obtained in high-resolution mode through the broadband 
filters f475 (not U Cam), f606, and f8 14 centred at 4760 A (width 
1458 A), 5907 A (width 2342 A), and 8333 A (width 2511 A), 
respectively, using the coronographic option. R Scl, where the 
detached shell is extended (diameter «40"), the large corono- 
graphic spot (3" in diameter) was used to effectively block the 
light of the bright star. For U Cam the angular size of the shell 
is smaller (diameter « 15") and the smaller coronographic spot 
(1'.'8 in diameter) was used. R Scl was observed in two rotation 
angles of the satellite (separated by 30°) in the f606 and f814 
filters to cover more of the shell. Finally, for each target star a 
nearby reference star with similar magnitude and colour was ob- 
served in the same mode (HD1760 for R Scl and HD20797 for 
U Cam). These images were used to subtract the point spread 
function (psf) from the images of the target sources. The expo- 
sure times are given in Table [2] In addition, short exposures of 
0.1s were obtained for the target stars without coronograph in 
all filters to provide estimates of the stellar fluxes. 



Table 1. Stellar properties [data from lSchoier et all d2005l) l. 



Source 


Var. type 


P 


D 




7* 






[days] 


[pc] 


[L©] 


[K] 


RScl 


SRb 


370 


290 


4300 


2630 


UCam 


SRb 


400 


430 


7000 


2700 


le 2. Observational parameters. 


Source 


Filter 




AA 


Exp. 








[A] 


[A] 


M 




RScl 


f475 


4760 


1458 


1734 






f606 


5907 


2342 


1206 a 






f814 


8333 


2511 


1206 a 




HD1760 


f475 






820 






f606 






330 






f814 






330 




UCam 


f606 






1060 






f814 






1060 




HD20797 


f606 






1050 






f814 






1050 





a Observed with this exposure time in two rotation angles of the satel- 
lite 



3. Data reduction 

We used the bias-corrected, flat-fielded, and cosmic-ray- 
corrected images supplied by the HST pipeline. For each image 
the background level was estimated and subtracted. The images 
were rebinned to a pixel size of 07026. At this point the rele- 
vant reference star images were subtracted from the target star 
images. This was achieved by shifting in position and multiply- 
ing the flux scale of the reference star image until the best result 
was obtained in a box (with a side of 13" and 6" for R Scl and 
U Cam, respectively) centred on the target star. Note that this 
subtraction procedure introduces some uncertainty in the radial 
behaviour of the scattered light distribution. 

The resulting (psf-subtracted) images were geometrically 
corrected and rotated to the (a, (5)-system using the correc- 
tion files on the ACS homepage and the relevant procedure in 
PYRAH3 (notably 'pydrizzle'). The images of R Scl obtained 
at the two rotation angles of the satellite (only for the f606 and 
f814 filters) were finally combined (the position of the target 
star was estimated in both images and they were appropriately 
shifted in position before combining). In this process the arti- 
facts produced by the l'/8 coronographic spot and the corono- 
graphic 'finger' were eliminated since these features appear in 
different regions of the rotated images. Finally, all images were 
median filtered. The resulting images are shown in Figs[TJand|2] 

The stellar fluxes were obtained from 0. 1 s exposure images. 
In the f475 filter the stellar flux density was obtained directly 
from the image. For the f606 and f814 filters this procedure was 
not possible, since the area on the CCD which contains most of 
the stellar flux was saturated even with 0.1 s of exposure time. 
The stellar fluxes therefore were estimated through psf fitting. 
The publicly available program 'Tiny Tim' was used to produce 
model HST psfs, which were rotated and resampled to have the 
same rotation and pixel scale as the observed images. This model 
psf image was fitted to the target image, excluding the central 
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saturated region, which made it possible to estimate the stellar 
flux densities in both filters. 

The presence of a circumstellar medium is apparent in 
all images. A comparison with the data obtained for R Scl 
bv IGonzalez Delgado et all d2001l) using EFOSC1 on the ESO 
3.6 m telescope is instructive. Clearly, the superior quality of 
the HST psf is crucial for this type of observations of extended 
low-surface-brightness emission in the vicinity of a bright ob- 
ject. Only an area with a diameter * 7" around the central star is 
seriously contaminated in the HST images. 



4. Results 

4.1. The nature of the circumstellar emission 

Circumstellar envelopes around AGB stars are cool objects (the 
kinetic and dust temperatures are <100K if we exclude the in- 
ner regions) and their emission, whether in lines or continuum, 
is very weak at optical wavelengths. Therefore, the circumstel- 
lar emission seen in our images must be due to scattered light. 
Further, it is scattered stellar light, as opposed to e.g. the case of 
the extreme carbon star IRC+10216, where the external parts of 
the CSE is see n through scattering of the interstellar radiation 
field (see e.g., iMauron &Hugginsl 119991). These conclusions 
were firmly drawn already by Gonzalez Delgado et al.l d2001l) in 
their analysis of the circumstellar optical emission detected by 
them towards the carbon stars R Scl and U Ant. 

The dominant scattering agent responsible for the emission 
in the images presented in this paper must be dust particles for 
the following reason. The width of a circumstellar line from a 
carbon star is about 30kms (a well establishe d resul t from 
radio line observations; e.g., ISchoier & Olofssonl d200ll) ). and 
this corresponds to a; 0.6 A at 6000 A. Consequently, a line is 
effectively diluted by a factor of about 4000 in the broadband 
filters used here. This does not absolutely exclude line scatter- 
ing. However, as will be shown below, for RScl the surface flux 
densities (i.e., flux per unit w avelength) in our images are co n- 
sistent with those obtained bv lGonzalez Delgado et al.l d2001l) in 
50 A wide filters centred on the resonance lines of Na and K. 
Specifically, the flux densities are the same in the f606 filter used 
by us and the 50 A wide f59 filter c entred on the NaD lines used 
by IGonzalez Delgado et al.l (1200 ll) . and therefore any contribu- 
tion from the NaD line to the f606 image must be negligible. 

Nevertheless, scattering in the reso nance lines of Nal and 
KI was observed in detached shells by iMaercker et alj d2010l) 
and in "normal" circum stellar envelopes of carbon stars by e.g. 
iGustafsson et al.l dl997l) . using spectroscopic techniques. 



4.2. The appearance of the circumstellar emission 
4.2.1. R Scl 

Due to the large extent of the circumstellar emitting region 
around R Scl compared to the field of view of the ACS in high- 
resolution mode only a part of it could be covered in the obser- 
vations, about 50% of the area and about 30% of the circum- 
ference. There is a marked limb-brightening, which follows an 
essentially circular arc, suggesting an overall spherically sym- 
metric ge ometry. This is consisten t with previo us observations 
of R Scl dGonzalez Delgado et al.L |2001L I2003I) and other de- 
tached shells. However, a circular arc fitted to the limb bright- 
ening does not have a centre at the stellar position (as is evident 
in Fig. [TJ, but it is offset by « ft' 7 at a position angle (PA) of 
about 160°. The same situation was found for the detached shell 



around TT C yg, where the shell ce ntre is offset l'/7 from the stel- 
lar position dOlofsson et al.L 120001) . This could be due to binary 
companions. We caution that for R Scl we cover only about 30% 
of the shell limb, and for this reason we use the stellar position 
as the centre of the shell in the following analysis. 

In order to analyse the circumstellar scattered light, az- 
imuthally averaged radial profiles (AARPs) were made in five 
PA intervals (see Fig. [TJ. The PA intervals are centred on 
355° (PA1), 32° (PA2), 57° (PA3), 77° (PA4), and 98° (PA5), 
and cover a range of +10°. Figure [3] shows the R Scl AARPs 
of the brightness distribution in PA4 in the f606 and f814 fil- 
ters, and PA5 in the f475 filter (PA4 is not entirely covered in 
the f475 filter). Also shown is the average AARP over the en- 
tire image for each filter. Overall, the surface brightness in the 
AARP covering the entire image is uniform out to a rather sharp 
outer edge. A low scattering optical depth means that the front 
and the rear parts of the CSE contribute equally to the result- 
ing brightness distribution (assuming isotropic scattering). The 
amount of scattered light is much lower in the f474 filter due 
to the marked drop in the light of carbon stars at shorter wave- 
lengths. The discussion here is therefore centred around the f606 
and f814 images. 

The limb brightening is apparent in the AARPs, and it is 
consistent with a CSE that has the form of a geometrically thin 
shell, and is in accordance w ith the polarimetric imaging of 
IGonzalez Delgado et al.l (120031) . which clearly shows the pres- 
ence of a detached shell. The limb brightening is much more 
prominent in the AARP of PA4 compared to the total AARP, i.e., 
it varies in strength with the PA, partly due to an actual strength 
variation but also because the size of the shell varies somewhat 
with the PA (mainly an effect of the shell not being centred on 
the star). There is weak emission extending at least 3" beyond 
the peak of the shell emission. The impression is that this is scat- 
tered light in material that protrudes in the radial direction from 
the shell, but the limited field of view of the ACS in high reso- 
lution mode combined with the large extent of the shell means 
that this area is poorly covered in the observations. There are in 
addition artifical radial structures, 'striatums' due to the corono- 
graphic mask, at low surface brightness. 

Looking in more detail at the brightness distributions reveals 
considerable patchiness, with a structure that is very similar in 
the f606 and f814 images. Indeed, in an image formed by the 
ratio of these two images there is no discernible patchiness, only 
noise around the average ratio. The same result is obtained when 
ratios involving the f475 image are produced. This also shows 
that there are very little colour variations over the shell in the 
scattered light, except that the f475/f606 and f606/f8 14 ratio im- 
ages show a systematic trend of the emission becoming redder 
with increasing distance from the star, most likely an effect of 
the wavelength dependence of the scattering asymmetry (see be- 
low). However, we caution that the reference star subtraction 
may affect any radial distribution. There appears to be no large- 
scale dominant pattern, but both the front and the rear part of the 
shell contribute to the scattered light, which complicates matter. 
The patchiness is only a weak modulation in the AARP aver- 
aged over all PAs, the peak-to-peak value of these variations be- 
ing less than +5%, while the patches affect the AARPs more at 
the separate PA intervals. The relevance of this for the detailed 
density distribution in the shell is discussed further in Sect. 14.: 



4.2.2. U Cam 

In this case, the ACS in high-resolution mode covers the full ex- 
tent of the circumstellar shell (Fig.|2j. However, the considerably 
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20 10 -10 20 10 -10 

offset [arcsec] offset [arcsec] 

Fig. 1. HST/ACS images of R Scl in the f475 (top left), f606 (top right), and f8 14 (bottom left) niters. The size of the coronographic 
spot used is 3". The l'/8 coronographic spot and the coronographic 'finger' is seen in the f475 image. About half of the shell is 
covered in these images. In the bottom right panel the sectors indicate the different PA intervals over which the AARPSs have been 
averaged. The centres of the sectors are all at the stellar position. 



smaller extent means that the artefacts caused by the central star 
light covers a larger fraction of the circumstellar emission. In ad- 
dition, since observations were obtained at only one rotation an- 
gle of the satellite, the artifact due to the 3" coronographic spot 
can be seen at the north-western edge of the CSE. Nevertheless, 
a limb brightening is clearly visible and the deviations from a cir- 
cle of this is less than 10%, i.e., it is consistent with a CSE that 
has the form of a geometrically thin shell with overall spherical 
symmetry. 

As for R Scl we create AARPs of the circumstellar emission 
around U Cam over the entire image and in different position 



angles. In this case the PA intervals had a width of 120° and are 
centred on 60° (PA1), 180° (PA2), and 300° (PA3) (see Fig.©. 
The region affected by the 3" coronographic spot is not included 
in PA3. Figure |4] shows the resulting AARPs for PA1 and the 
AARPs over the entire image in the f606 and f814 filters. 

Similarly to the case of R Scl, the surface brightness in the 
AARP of the entire image is uniform out to a rather sharp outer 
edge. The limb brightening, however, is much more apparent 
here (particularly in the f606 image). 

There is considerable patchiness in the brightness distribu- 
tions, which is fully consistent in the f606 and f8 14 images. 
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Fig. 2. HST/ACS images of U Cam in the f606 (top left), and f814 (top right) niters, and the different PA intervals over which the 
brightness distributions were averaged (bottom left). The size of the coronographic spot used is 1"8. The 3" coronographic spot is 
seen at PA « -45° at the edge of the shell. 



There appears to be no large-scale dominant features, but the 
outer edge clearly exhibits a kind of undulating pattern, which 
may simply be due to a random distribution of varying surface 
brightness in the shell. 

4.3. Scattering in detached dust shells 

iMaercker et al.l d2010l) analysed the scattered light around the 
carbon star U Ant by fitting calculated radial intensity distribu- 
tions to the observed AARPs using the optically thin approxima- 
tion for scattering of stellar light by dust grains in geometrically 
thin circumstellar shells. They assumed that the grains are of a 



single size (0.1 yum), consist of amorphous carbon, and ha ve a 
density of 2 g cm 3 . The optical constants were taken from lSuhl 
(l2000h - who analysed the spectral energy distributions (SEDs) 
of carbon stars. The scattering efficiency (as a function of scat- 
tering angle) was calculated with the Mie theory. The radial den- 
sity distribution of the dust wass assumed to follow a Gaussian 
distribution centred at a radius R with a full wi dth at half max- 
imum AR. We follow the same procedure as in IMaercker et al.l 
(1201 Oh to derive the size and width of the detached shells around 
R Scl and U Cam. Calculated brightness distributions are fitted 
to the AARPs in the different PA intervals for R Scl and U Cam. 
Examples of the fits are shown in Figs. [3] and [4] Here we adopt 
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Distance to the star [arcsec] 




Distance to the star [arcsec] 




5 10 15 20 

Distance to the star [arcsec] 



Fig. 3. HST/ACS AARPs of R Scl in the f475 (left), f606 (middle), and f8 14 (right) filters. The panels show the AARP in PA4 (solid 
line) (for the f475 filter PA5 is shown), the AARPs averaged over all PAs (dotted line), and fits assuming scattering of stellar light 
in a detached shell of dust (dashed line). See text for details. 



Table 3. Results of the fits to the AARPs in different PA intervals 
for R Scl and U Cam. 



filter 


PA 


R 


AR 


^peak 




interval 


["] 


["] 


[ergs" 1 cm" 2 A -1 "~ 2 ] 


RScl 










f475 


PA2 


18.5 


0.9 


0.7 x 10~ 17 




PA5 


19.0 


0.9 


0.4 x 10- 17 


f606 


PA1 


18.7 


1.8 


6.3 xl0~ 17 




PA2 


18.9 


1.6 


6.2 x 10~ 17 




PA3 


19.5 


1.1 


5.5 xlO- 17 




PA4 


20.0 


0.8 


5.0 x 10~ 17 




PA5 


19.3 


1.6 


4.8 x 10~ 17 


f814 


PA1 


18.8 


1.8 


6.0 x 10 17 




PA2 


18.9 


1.6 


5.9 x 10~ 17 




PA3 


19.6 


1.1 


5.1 x 10~ 17 




PA4 


19.9 


0.7 


4.8 x 10- 17 




PA5 


19.4 


1.9 


4.6 x 10~ 17 


UCam 










f606 


PA1 


7.6 


0.5 


1.3 x 10~ 17 




PA2 


7.7 


0.4 


1.2 x 10~ 17 




PA3 


7.7 


0.7 


0.9 x 10- 17 


f814 


PA1 


7.6 


0.6 


3.6 x 10~ 17 




PA2 


7.7 


0.5 


2.9 x 10~ 17 




PA3 


7.6 


0.8 


2.8 x 10- 17 



a grain size of 0.15 /mi, because the higher forward scattering 
efficiency of these grains gives a better fit to the AARPs for both 
R Scl and U Cam (see below). 

In general, the peak due to limb brightening and the tail could 
be fitted reasonably well, giving good estimates of the shell ra- 
dius and width. The results for R and AR at the different PA 
intervals for R Scl and U Cam are given in Table [3] On average, 
the dust shell radii for R Scl and U Cam are 19'.'2 (correspond- 
ing to a linear radius of 8x10 cm) and 7'.'7 (corresponding to 
a linear radius of 5xl0 16 cm), respectively. The average dust 
shell widths based on the fits are 1'.'2 (corresponding to a linear 
width of 5xl0 15 cm) and ft' 6 (corresponding to a linear width 
of 4xl0 15 cm) for R Scl and U Cam, respectively. Modelling 
of the polarised scattered flux from R Scl gives a shell r adius 
of 20", and a width of 2" (iGonzalez Delgado et alll2003l) . The 
ground-based observations are subject to seeing, which has the 
effect of moving the peak to larger radii and widening the shell, 



and so the results are consistent with each other. ISchoier et alJ 
J2005I) derive a dust shell radius of «27" by modelling the 
SED of R Scl, but with considerable uncertainty («+14"). A 
recent SMA CO(/ = 2-l) interferometer map of the R Scl shell 
gives an estimated CO shell radius of « 18" at a resolution of 
about 4" (Dinh-van-Trung, priv. comm.). For U Cam the size 
of the CO shell determined in IRAM PdB interferometer maps 
is 7" 3, which excellen tly agrees with the results derived here 
dLindqvist et al.L fl999t) . The width is more difficult to derive in 
the interferometer data, but is estimated to be =1". Thus, we 
conclude that in both cases the dust and gas shells coincide spa- 
tially within the uncertainties (< 1" and « 2" in the case of R Scl 
and U Cam, respectively). This contrasts with the finding for the 
detached shell(s) around the carbon st ar U Ant, where the gas 
and dust shells have clearly separated ( Gonzale z Delgado et all 
l2003blMaercker et all 12010) . 

The peak brightness of the limb brightening is given for the 
different PA intervals in Table[3] It is clear that the limb brighten- 
ing is much less pronounced in the AARP averaged over all PAs 
for both R Scl and U Cam. For R Scl the fits to the AARPs in the 
different PA intervals show that the estimated radius of the shell 
varies by « +1" (caused by the centre of the shell most likely not 
being centred at R Scl). At the same time the width of the shell 
is * 1"2, with the consequence that the limb brightening in the 
AARP over the entire image is smeared out. The estimated ra- 
dius of the shell around U Cam is essentially the same in all PA 
intervals, but the width varies considerably, which also leads to 
a smearing of the limb brightening in the AARP averaged over 
all PAs. 

The scattering of the grains occurs preferentially in the for- 
ward direction when the grain size becomes comparable to the 
observed wavelength [the deta ils of the scatter i ng an d its angu- 
lar dependence is described in iMaercker et all (1201 Oh l. and this 
explains that the scattered light does not become significantly 
weaker closer to the star which would otherwise be expected 
for isotropic scattering in a geometrically thin shell. Indeed, this 
property can be used to estimate the size of the grains because 
their forward scattering efficiency depends on this (as is shown 
in Fig. [5] through g, which is a measure of the directional depen- 
dence of the scattering). However, it is clear that the simplified 
model does not fully explain the AARPs at line-of-sights closer 
to the star. We find a reasonable fit to the AARPs for a grain size 
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Table 4. Quantitative results for R Scl and U Cam from the fits to the observed brightness distributions averaged over all PAs ". 



Source 


filter 


R 


AR 


F, 


^CSE 


^CSE 


S CSE 


Fcse/F, 








["] 


["] 


[erg s cm -2 ] 


[erg s~' cirT 2 A -1 "~ 2 ] 


[erg s" 1 cm -2 ] 


[erg s cm -2 A -1 ] 






RScl 


f475 


18.8 


0.9 


4.9 x lO - * 


0.5 x 10 _1/ 


0.1 x 10-' u 


0.6 x 10~ 14 


17.0 x lO- 


4 




f606 


19.3 


1.4 


1.3 x 10~ 7 


4.9x 10~ 17 


1.4 x lO" 10 


6.0 x 10- 14 


ll. 1 x 10" 


4 




f814 


19.3 


1.4 


1.3 x 10- 7 


4.4 x 10~ 17 


1.4 x 10" 10 


5.5 x 10- 14 


10.3x10" 


4 


UCam 


f606 


7.7 


0.5 


1.4 x 10~ 8 


0.6 x 10~ 17 


0.4 x 10-" 


0.2 x 10~ 14 


2.6 x 10" 


4 




f814 


7.6 


0.6 


5.4 x 10" 8 


3.0x 10~ 17 


1.4 x 10" 11 


0.6 xlO" 14 


2.5x10" 


4 



The columns give the radius (R) and the FWHM (AR) of the Gaussian dust density distribution, the stellar flux (F,), the average surface 
brightness (ZcseX the A ux of the scattered stellar light (F C seX the A ux density of the scattered stellar light (ScseX an d the ratio between the 
circumstellar and stellar fluxes (Fcse/F,). See text for details. 
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U Cam - f606 
PA1 

'/ i\ 
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Distance to the star [arcsec] 
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5 6 7 8 9 
Distance to the star [arcsec] 

Fig. 4. HST/ACS AARPs of U Cam in the f606 (top) and f814 
(bottom) filters in the PA1 interval, the AARPs averaged over all 
PAs (dotted line), and fits assuming scattering of stellar light in 
a detached shell of dust (dashed line). See text for details. 



of 0. 15 pm, but there is a notable systematic underestimate in the 
f814 filter. A distribution of grain sizes would probably lead to 
a better fit in all filters. Also, the details of the AARPs are not 
well fit. This is likely due to the patchiness observed in the shells 
(i.e., the detailed density distribution of the dust). The average 
surface brightnesses of the shells in the different filters are there- 
fore obtained as straight averages over the area outside a radius 
of 5" and inside a radius of w 17" and 7" for R Scl and U Cam, 
respectively. The resulting values are given in Tables [3] (for the 



different PA intervals) and B ( for the entire image). For R Scl, 
IGonzalez Delgado et alj J2001I) obtained 2 x 1(T 17 and4x 1(T 17 
(in units of erg s _1 cm" 



-1 n-2. 



with some considerable uncer- 



tainty) in 50A wide filters centred at 590 nm and 770 nm, respec- 
tively (i.e., the wavelength ranges of these filters lie inside those 
of the f606 and f8 14 filters, respectively). Thus, within the un- 
certainty, the flux densities are the same in the narrow and the 
broad filters, suggesting that even in the narrow filters the con- 
tribution from atomic scattering (by Na and K in the 590 nm and 
770 nm filter, respectively) is limited, and the broad filters used 
on the HST are completely dominated by dust scattering. 

Of a more quantitative importance is the ratio between the in- 
tegrated surface brightness and the stellar flux, i.e., the fraction 
of the stellar light that is scattered. The integrated scattered light 
fluxes are obtained by summing the AARPs covering the entire 
image over all angles. The estimated uncertainty in the stellar 
fluxes due to the psf fitting is estimated to be a factor of as 2. The 
circumstellar and stellar fluxes as well as the resulting circum- 
stellar to stellar flux ratios are given in Table|4] The results show 
that on average the scattering process is very optically thin. 

4.4. Dust shell masses 

We proceed with a simple analysis of the observed scattered 
emission based on the assumption of optically thin dust scat- 
tering in a geometrica l ly thi n shell following the work of 
IGonzalez Delgado et al.l (1200 ll) . The circumstellar scattered flux 
(at the distance of the object), integrated over the total scattering 
wavelength range (defined by filter f) in the optically thin regime 
is given by 



F sc ,t = - — N sc f ct sc ,aF*a 



cU. 



(1) 



where F t ,\ is the stellar surface flux. Hence, the expected circum- 
stellar to stellar flux ratio is given by 



F*, f 4ttR 2 



1 I °~ sc,aF t/ id/l 

-Nj f 



2 sc 
sh 



(2) 



where 7V SC is the number of scatterers, <x sc the scattering cross 
section, and R s h the shell radius. 

For dust scattering, and assuming single-sized, spherical 
grains, we have the following relations for the number of scat- 
terers and the scattering cross section 



/Vsc 



3M d 
4np i a i 



O" SC./1 = Qsc.A^a , 



(3) 
(4) 



8 



H. Olofsson et al.: High-resolution HST/ACS images of detached shells around carbon stars 



where M& is the dust mass in the shell, p g the density of a dust 
grain, a the grain radius, and gsc,,) the grain scattering efficiency 
(assuming isotropic scattering). 2sc,i i s calculated with the Mie 
theory and the optical constants from ISuhl d2000h for carbona- 
ceous dust grains with a radius of 0.15 fim and a grain density of 
2gcrrT 3 . 

Qsc.a varies strongly with wavelength and grain size, see 
Fig. [5] Consequently, Q sc ,a varies considerably over the wave- 
length range covered by one filter. A directional dependence 
of the scattering, which depends on the grain size, further 
complicates the problem. As opposed to the analysis of the 
AARPs we here assume isotropic scattering for reasons of 
simplicity. Finally, the assumption of one constant grain size 
throughout the shell adds an additional uncertainty. This re- 
sults in significant uncertainties in the determined dust masses 
in the shell. Together with the uncertainty in the stellar fluxes 
(and hence Fqse/F, ratios), we estimate the resulting dust 
masses to be uncertain by a factor »5. The resulting shell 
masses based on the measured circumstellar to stellar fluxes 



for R Scl and U Cam are 3xl(T 6 M and 3xl(T 7 M , 



respec- 



tively (Table Hi. For R Scl. iGonzalez Delgado etail (120031) de- 
rived a dust mass, ba sed on models of pol arised scattered light, 
of 2xl(T 6 M , while ISchoier et alJ (120051) derived a dust mass 
of (3.2 + 2.0)xl( T 5 Mo by modelling the dust thermal emis 



sion. For U Cam ISchoier et al.l ( 20051) deri ved only an upper 
limit to the dust mass. Sch oier et al T (l2005h derived gas shell 
masses, based on CO radio line modelling, of 2.5xl0~ 3 M and 
lxl0~ 3 M for R Scl and U Cam, respectively, leading to dust- 
to-gas mass ratios of lxlO 3 and 3 xl0~ 4 , rather low values for 
carbon stars dRamstedt et al.Ll2008h . 

4.5. The small-scale structure of the circumstellar medium 

The detached, geometrically thin shells observed here are ex- 
cellent for studies of the detailed density distribution of the cir- 
cumstellar medium, in particular in dust-scattered stellar light, 
where the optical depth is low. For both R Scl and U Cam sig- 
nificant patchiness is observed in the images, and it is close to 
identical in the different filters for the respective object. Due to 
the smaller size of the shell the patchiness is less obvious for of 
U Cam. In particular the artefacts due to the psf-subtraction in 
the inner parts of the shell make it difficult to analyse the U Cam 
data. 

In order to analyse the small-scale structure in more detail, 
unsharp masks (smoothing with a Gaussian having a FWHM 
of 3.5") were subtracted from the f814 images for R Scl and 
U Cam, Fig. [6] In these images the "dumpiness" is much more 
apparent. Radially averaged azimuthal profiles (RAAPs) were 
created from the sharp images at three radial distances for R Scl 
(1 1'.'5, 15'.'5, and 1971) and two radial distances for U Cam (6'.'1 
and 7'.' 4). The azimuthal ranges covered by the RAAPs for R Scl 
and U Cam are indicated in Figs. [JJ and [2] respectively. The 
RAAPs were averaged over a radial range corresponding to 5% 
of the radial distance. The resolution along the azimuthal direc- 
tion is 07025 for R Scl and 072 for U Cam. The resulting RAAPs 
are shown in the upper panels of Figs. [7] and [8] for R Scl and 
U Cam, respectively. The patchy structure can clearly be seen in 
the RAAPs of R Scl, and the profiles at the different radii show 
the same general type of structure, indicating that the clumps 
are evenly distributed over the shell. The U Cam image is nois- 
ier than the R Scl image, and the artifacts from the coronograph 
cover a large fraction of the shell, but the patchy structure can 
clearly be seen also here in the RAAP at 7'.'4. The lower panels 
show the corresponding Fourier transforms of the RAAPs. They 




0.4 0.6 0.8 

Wavelength [jam] 



Fig. 5. The scattering efficiency Q sc ,i (thick lines) and a mea- 
sure of the directional dependence of the scattering g (thin lines) 
as functions of wavelength for 0.1 /i, 0.15 jum, and 0.2 //m sized 
grains (solid, dashed, and dotted lines, respectively). 

Table 5. The derived dust shell-masses based on the observed 
Fcse/F* ratios. 



Source 


filter 


M d 
[1(T 6 M ] 


M 5 " 
[10- r M o ] 


RScl 


f475 


2.6 


2.5 




f606 


2.1 






f814 


4.0 




UCam 


f606 


0.2 


0.9 




f814 


0.3 





a The shell gas mas s estimated from CO radio line modelling 
( ISchoier e7aUl2005h . 



show that some tendencies for structures arise on scales > 1" in 
both cases. 

The patches seen in the images are here interpreted as signs 
of a small-scale structure in the form of clumps in the detached 
shells. We made a crude estimate of the "dumpiness" in the 
R Scl shell through randomly distributing clumps of an equal 
size in a shell of a radius 1972 and a width 172. Each clump 
produces optically thin scattering with a directional dependence 
according to the grain properties used in Sect. 14.31 A reasonable 
resemblance with the observed images and AARPs is found for 
2000 clumps of a size 079 (the FWHM of a Gaussian clump). 
Structures on scales > 1" arise because of overlapping clumps. 
As an example, the simulated image in the f814 filter is shown 
in Fig. [9] The RAAPs, in the same radial and azimuthal ranges 
used for Fig. [7] obtained from the sharpened version of Fig. [9] 
(using the same procedure as for the observed image) are shown 
in Fig [TO] together with the associated frequency distribution. 
These are also similar to the observed ones except for a lack of 
structure on s maller scales, probably an effect of a distr i bution 
of clump sizes. [Bergman et al.1 d 1993b and lOlofsson et alJ d2000h 
estimated that the detached CO shells of the carbon stars S Set 
and TT Cyg are composed of about 1000 clumps of a size of 
about 075 (FWHM), in line with the findings here. Adopting a 
value of 2000 clumps leads to an estimated clump dust mass of 
about 2xl0- 9 M o for RScl. 
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Fig. 6. Images in the f8 14 filter of R Scl (right) and U Cam (left) after the subtraction of an unsharp mask, emphasising the structure 
in the detached shells. 
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Fig. 7. Top panel: RAAPs at three distances from R Scl: 1971 
(top), 15'.'5 (middle), and 1175 (bottom). The azimuthal profiles 
are averaged over a radial range of 5% of the radius. Bottom 
panel: Corresponding power spectra for the different profiles (in 
the same order from top to bottom as in the upper panel). The 
profiles in both panels are offset for clarity. The original zero 
levels are indicated by the dashed lines. 




frequency [arcsec 



Fig. 8. Top panel: RAAPs at two distances from U Cam: 774 
(top) and 67 1 (bottom). The azimuthal profiles are averaged over 
a radial range of 5% of the radius. Bottom panel: Corresponding 
power spectra for the different profiles (in the same order from 
top to bottom as in the upper panel). The profiles in both panels 
are offset for clarity. The original zero levels are indicated by the 
dashed lines. 



A comparison of Figs Q] [2] and [6] with Fig. [9] shows a dif- 
ference between the simulations and the observations. The ob- 



servations show sections of possible double-shell structures or 
features, which appear to be fragmented shells or possibly fila- 
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merits, as opposed to the simulations where the clumps are ran- 
domly distributed in the shell. It is unclear if these structures are 
the results of wind-wind interactions due to instabilities, or of 
density or velocity inhomogeneities in any of the two winds, or 
both. 



5. Discussion and conclusions 

5.1. Characteristics of the detached shells 

At the distances of R Scl and U Cam the angular radii of the 
shells corr espond to linear radii of 8 x 1 1 6 cm and 5xl0 16 cm, re- 
spectively. ISchoier et alJ (120051) estimate that the expansion ve- 
locities for the gas shells of R Scl and U Cam are 15.5 km s~' 
and 23.0kms~ 1 , respectively, and, since there appears to be no 
separation between the gas and the dust shells, the correspond- 
ing dynamical ages of the shells are consequently about 1700 yr 
(R Scl) and 700 yr (U Cam). Thus, these are the youngest de- 
tached shells observed in detail so far. The average shell widths 
convert into linear widths of 5xl0 15 cm and 4xl0 15 cm for R Scl 
and U Cam, respectively, and the corresponding time scales are 
100 yr (R Scl) and 50 yr (U Cam). The shell width/radius ra- 
tios lie in the range 0.05 - 0.1. The estimated dust and gas shell 
masses are 3xlO~ 6 M and 2.5xl0~ 3 M Q , respectively, forR Scl 
and 3xl0~ 7 M o and lxlO~ 3 M , re spectively, for U Cam . To 
this we can add the recent result by iMaercker et al.l d2010t) for 
the carbon star U Ant, where the shell has a dynamical age of 
2700 yr, a width corresponding to 140 yr, and the dust and gas 
masses are estimated to be 5xl0 _5 M o and 2.5x l0~ 3 M H . Our 
new d ust mass estimates agree with the finding by ISchoier et alJ 
(120051) that the shell masses increase with the age of the shell, 
which is possibly an effect of interaction between the shell and 
a previous slower stellar wind. The dust-to-gas ratio also in- 
creases with the age of the shell (3xl0~ 4 , lxl0~ 3 , and 2xl0" 2 
for U Cam, R Scl, and U Ant, respectively). This is curious, be- 
cause on e expects the dust to g radually leave the gas shell as it 
expands (IMaercker et all 1201 Oh . at least if the wind is smooth, 
dumpiness, possibly arising as the shell expands, and/or wind- 
wind interaction may be an explanation here as this will affect 
the gas-grain drift and may lead to accumulation of material. 



5.2. Origin of the detached shells 

It has been argued by various authors that the geomet- 
rically thin, detached shells found around carbon stars 
are the effects of mass-loss rate modulations during ther- 



mal pulses (Olofsson et al., 1990, 


19961 ISchroeder et al. 


19981 Gonza 


ez Delgado et all 12003 


: ISchoier et al. . 2005: 


Mattsson et al. 


. 120071 Maercker et all 1 


2010). In particular, the 


detailed models of Mattsson et alJ (120071) show that aeometri- 



cally thin shells appear under certain circumstances as an effect 
of a mass-loss eruption and a subsequent interaction with a pre- 
vious slower wind. The results on R Scl and U Cam presented 
here are consistent with such a scenario. They add further sup- 
port to the presence of highly isotropic mass loss during the shell 
formation event, and, if interaction plays a major role, this must 
apply also to the mass loss prior to the eruption. If this interpre- 
tation is correct, both stars must be in the aftermath of a very re- 
cent thermal pulse (<2000 years ago). The shell formation time 
scale is as short as < 100 yr. This is an estimate that is likely 
not affected by interaction with a previous slow wind because of 
the relative youth of the shells. The amount of mass ejected is 
relatively small, w 10~ 3 M Q . 
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Fig. 9. A simulated image of the detached shell around R Scl 
seen in dust-scattered light in the f814 filter (see the text for de- 
tails). 



An estimate of the mass-loss rate during the formation of 
the shells is obtain ed using the shell gas mass estimates of 
ISchoier et all (120051) and the time scales estimated here. The re- 
sults are 3xlO~ 5 and 2xl0~ 5 M Q yr 1 for R Scl and U Cam, re- 
spectively. This ignores any effects of e.g sweeping up of mate- 
rial. An estimate of the mass-loss rate prior to the shell ejection 
is difficult to give. We limit ourselves here to an estimate of the 
density contrast between the shell and the surrounding medium. 
In the light of optically thin scattering, the density contrast be- 
tween the shell and the medium outside the shell is at least a 
factor of 30 in the case of R Scl (for which the data has the high- 
est S/N). This should be a rough estimate of the mass-loss rate 
contrast as well. The mass-loss rate following the ejection of the 
shell is difficult to estimate from our data, because it depends 
critically on the psf subtraction and the propertie s of the grain 
scatter ing (in particular its angular dependence). ISchoier et al.l 
(I2005h estimated, based on CO radio line modelling, that the 
present mass-loss rates are 3xl0~ 7 and 2xl0~ 7 M yr~' for R Scl 
and U Cam, respectively. In conclusion, the mass-loss rate mod- 
ulations are substantial, of the order of two magnitudes. 

5.3. Origin of clumps in the detached shells 

The small-scale structure of the circumstellar medium is of in- 
terest for various reasons, e.g., it will affect the radiative trans- 
fer and hence the estimates (like those of mass-loss rates and 
molecular abundances) obtained from modelling molecular line 
emission, it will affect the circumstellar chemistry (incl. the pho- 
todissociation), and it will affect the efficiency of the mass loss. 

The observational evidence for dumpiness in AGB CSEs 
is scattered and not very conclusive. High angular resolution 
near-IR observations certainly suggest a very inhomogeneous 
medium close to the st ar, where even proper motions of clumps 
have been measured (IWeigelt et all 120021) . Far-infrared and 
radio observations of the circumstellar medium normally do 
not have the resolution to identify any small-scale structure. 
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Fig. 10. Same as Fig. [7] but with the data input from the sharp- 
ened version of the simulated image in Fig. [9] (see text for de- 
tails). 



However, observations of in particular SiO and H2O maser emis- 
sion suggest the presence of clumps also in the CSEs. The very 
long baseline interferometry (VLBI) observations of SiO masers 
towards AGB stars suggest s pot sizes as small as a few 10 12 cm 
in regions close to the star (IColomer et al 1 119921) . Indeed, for 
reasonable mass-loss rates the gas within 10 14 cm of the star 
must be substantially clumped to provide the densities required 
for the SiO masers (about 10 cm -3 ). Also H?Q masers show 
spot sizes as small as a few 10 cm dlmai etall fl997b . The 



characteristic OH maser spot size appears to be a few 10 cm 
at a radius of a few 1 15 cm, i.e., further out in the envelope 
dChapman et all |1994|) . Note though that maser observations 
are difficult to interpret in terms of a detailed density structure, 
and they tend to enhance any density contrasts present. Larger- 
scale patchiness in CO radio line brightness maps of detached 
shells have been int erpreted as a consequence of a large num- 
ber of small clumps dBergman et all[i993l : l01ofsson et al.Lll996l 
[2000h . 

On the other hand, high resolution optical observations of 
planetary nebulae (PNe) show conclusive evidence of small- 
scale structure, in particular, the cometary globules seen in 
the Helix nebula. These are about 1" in size and are clearly 
affecte d by the strong radia t ion field from the central whit e 
dwarf dMeaburn et all ,119961 |1998t lO'dell & Handronl 11996b . 



iHuggins et al.l d 19921 l2002h detected CO radio line emis- 
sion from these clumps and estimate their masses to be 
about 10~ 5 Mq. It is likely that s i milar st r uctures are present 
also i n other PNe dWilsonl 119501: IZanstraL Il955t ISpeck et all 
20031): the Helix happens to be one of the most nearby PNe. 



Hug gins & Mauronl d2002l) studied the small-scale structure in 



two well-known CSEs, those around the carbon star IRC+10216 



and the young PN NGC 7027, with the aim to identify at which 
point the small-scale structure arises. They found no evidence 
of dumpiness in any of the two envelopes. We note in pass- 
ing that clumps have also been observed in expanding shells 
aroun d stars in other circumstances . Examples are O ass ocia- 
tions dFriemanl 1 19541) . HII regions dReipurth et all 1 1997b, and 
in the environments of y oung stellar objects ( Rudolph & Welchl 
ll988tlKeto&Holll989t) . 

The combination of optically thin scattering due to dust and 
the narrow widths of the detached shells is excellent for study- 
ing the small-scale structure. A potential drawback is that the 
mass-dominant component, the ga s, may not be d i stribu ted in 
the same way as the dust (see e.g.. iMaercker et all 1201 Oh . The 
R Scl sharpened image lends itself to a quantitative analysis of 
the clump properties. It suggests clump sizes of about 0/9 in di- 
ameter, corresponding to linear clump diameters of 4xl0 15 cm. 
There is an estimated total of about 2000 clumps in the shell. The 
average clump dust mass is consequently about 2xl0~ 9 M Q . The 
estimated dust-to-gas ratio suggests a clump gas mass of about 
2xlO~ 6 M Q . This is somewhat lower than the mass estimat es of 
the Helix nebula cometary globules dHuggins et all 120021) . but 
the uncertainties in both estimates are substantial. 

A likely explanation for the creation of detached shells 
around carbon stars is an increased mass-loss during the He- 
shell flash, with a subsequent int eraction of a faster wind run- 
ning into a previous , slower wind dSteffen & Schonberneit 120001 : 
ISchoier et all 120051: iMattsson et all 2007b . The question arises 
whether the clum py structure is alread y present in the ejected gas 
[as argued e.g. bv lDvson et al.l dl989t) l. or whether it emerges in 
the ex panding gas as a consequence of instabilities [as argued 
e.g. bv lCapriottil d 19731) for clumps in PNe]. 

It is difficult to judge whether the clumps reflect inhomo- 
geneities already present early, e.g. in the wind-accelerating 
zone at a distance of a few stellar radii from the star, or even 
in the stellar photosphere or below, as convective cells. We can 
identify two problems in this context. The clumps ejected from 
the star in this scenario are likely to have a distribution of prop- 
erties (mass, size, etc.). If so, one would expect clumps with dif- 
ferent properties to reach somewhat different terminal velocities, 
and once they start to interact with the previous slower wind 
(which is presumably also clumpy), one would expect clumps 
of different sizes or masses to accrete different amounts of ma- 
terial from the slow wind, and thus aquire different speeds for 
the second time. Thus it is hard to explain why the clumps 
seem to be located at equal distances from the star in a thin 
shell. On the other hand, as a result of the faster wind collid- 
ing with a slower wind, instabilities of various types may occur, 
e.g., Rayleigh-Taylor, Kelvin-Helmholtz, thermal, and gravita- 
tional instabilities. It is beyond the scope of this paper to anal- 
yse this in detail, but we note that numerical models indicate 
that structures roughly corresponding to the observed sizes and 
densities may indeed oc cur due to Rayleigh-Taylor instabilities 
dMvasnikov et all l2000). Once the clump characteristics are de- 
termined and their formation mechanism understood, the clump 
properties can possibly be used as a diagnostic of e.g. stellar 
evolution. 
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